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Abstract 

We construct a DGP inspired braneworld scenario where a scalar field non- 
minimally coupled to the induced Ricci curvature is present on the brane. First 
we investigate the status of gravitational potential with non-minimal coupling and 
observational constraints on this non-minimal model. Then we further deepen the 
idea of embedding of FRW cosmology in this non-minimal setup. Cosmological 
implications of this scenario are examined with details and the quintessence and 
late-time expansion of the universe within this framework are examined. Some 
observational constraints imposed on this non-minimal scenario are studied and 
relation of this model with dark radiation formalism is determined with details. 
PACS: 04.50.+h, 98.80.-k 
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1 Introduction 



Based on light-curves analysis of several hundreds type la supernovae [1,2], observations 
of the cosmic microwave background radiation by the WMAP satellite [3] and other CMB- 
based experiments [4,5], it has been revealed that our universe is currently in a period 
of accelerated expansion. Some authors have attributed this late-time expansion of the 
universe to an energy component referred to as dark energy. The simplest example in this 
regard is the cosmological constant itself which provides a model of dark energy. However, 
it is unfavorable since it requires a huge amount of fine-tuning [6]. Phantom fields [7], 
quintessence [8] and modification of gravitational theory itself [9,10] are other attempts to 
explain this late time expansion of the universe. In the spirit of the modified gravitational 
theory, Carroll et al have proposed R^ 1 modification of the usual Einstein-Hilbert action 
[11]. It was then shown that this term could give rise to accelerating solutions of the field 
equations without dark energy. 

On the other hand, theories of extra spatial dimensions, in which the observed universe 
is realized as a brane embedded in a higher dimensional spacetime, have attracted a lot 
of attention in the last few years. In this framework, ordinary matters are trapped on the 
brane but gravitation propagates through the entire spacetime [9,12,13]. The cosmological 
evolution on the brane is given by an effective Friedmann equation that incorporates the 
effects of the bulk in a non-trivial manner [14]. From a cosmological view point, the 
importance of brane models lies, among other things, in the fact that they can provide 
an alternative scenario to explain the late-time accelerating expansion of the universe. 

Theories with extra dimensions usually yield correct Newtonian limit at large dis- 
tances since the gravitational field is quenched on sub-millimeter transverse scales. This 
quenching appears either due to finite extension of the transverse dimensions [12,15] or 
due to sub-millimeter transverse curvature scales induced by negative cosmological con- 
stant [13,16-19]. A common feature of these type of models is that they predict deviations 
from the usual 4-dimensional gravity at short distances. The model proposed by Dvali, 
Gabadadze and Porrati (DGP) [9] is different in this regard since it predicts deviations 
from the standard 4-dimensional gravity even over large distances. In this scenario, the 
transition between four and higher-dimensional gravitational potentials arises due to the 
presence of both the brane and bulk Einstein terms in the action. In this framework, ex- 
istence of a higher dimensional embedding space allows for the existence of bulk or brane 
matter which can certainly influence the cosmological evolution on the brane. Even if 
there is no 4-dimensional Einstein-Hilbert term in the classical theory, such a term should 
be induced by loop-corrections from matter fields [20,21]. Generally one can consider the 
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effect of an induced gravity term as a quantum correction in any brane-world scenario. 
A particular form of bulk or brane matter is a scalar field. Scalar fields play an important 
role both in models of the early universe and late-time acceleration. These scalar fields 
provide a simple dynamical model for matter fields in a brane-world model. In the context 
of induced gravity corrections, it is then natural to consider a non-minimal coupling of the 
scalar field to the intrinsic (Ricci) curvature on the brane that is a function of the field. 
The resulting theory can be thought of as a generalization of Brans-Dicke type scalar- 
tensor gravity in a brane-world context. In contrast to common belief, the introduction 
of non-minimal coupling (NMC) is not a matter of taste; NMC is instead forced upon 
us in many situations of physical and cosmo logical interest. There are many compelling 
reasons to include an explicit non-minimal coupling in the action. For instance, NMC 
arises at the quantum level when quantum corrections to the scalar field theory are con- 
sidered. Even if for the classical, unperturbed theory this NMC vanishes, it is necessary 
for the renormalizability of the scalar field theory in curved space. In most theories used 
to describe inflationary scenarios, it turns out that a non- vanishing value of the coupling 
constant cannot be avoided. In general relativity, and in all other metric theories of grav- 
ity in which the scalar field is not part of the gravitational sector, the coupling constant 
necessarily assumes the value of | . The study of the asymptotically free theories in an 
external gravitational field with a Gauss-Bonnet term shows a scale dependent coupling 
parameter. Asymptotically free grand unified theories have a non-minimal coupling de- 
pending on a renormalization group parameter that converges to the value of | or to any 
other initial conditions depending on the gauge group and on the matter content of the 
theory. An exact renormalization group study of the A(/> 4 theory shows that NMC= | is 
a stable infrared fixed point. Also in the large N limit of the Nambu-Jona-Lasinio model, 
we have NMC= |. In the O(N)- symmetric model with V = A0 4 , NMC is generally 
nonzero and depends on the coupling constants of the individual bosonic components. 
Higgs fields in the standard model have NMC= or |. Only a few investigations produce 
zero value (for a more complete discussion of these issues we refer to papers by V. Faraoni, 
specially Ref. [22] and references therein). In view of the above results, it is then natural 
to incorporate an explicit NMC between scalar field and Ricci scalar in the inflationary 
paradigm and in quintessence models. In particular it is interesting to see the effect of 
this NMC in a DGP-inspired braneworld cosmology. 

There are several studies focusing on the braneworld models with scalar field [22-31]. 
Some of these studies are concentrated on the bulk scalar field minimally [23-25] or non- 
minimally [26-28] coupled to the bulk Ricci scalar. Some other authors have studied the 
minimally [29,30] or non-minimally [22,31] coupled scalar field to the induced Ricci scalar 
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on the brane. However, none of these studies have considered consequences of embedding 
of FRW cosmology with non-minimally coupled brane-scalar field into DGP scenario ex- 
plicitly. Specifically, none of these studies obtained five dimensional metric components 
and full dynamics of the braneworld within this scenario trivially. In addition, there are 
very limited literature on late time behavior and quintessence with non-minimally coupled 
scalar field and even these studies have not performed explicit and detailed calculation of 
the late-time dynamics. 

In this paper, in the spirit of DGP inspired gravity, we study the effect of an induced 
gravity term which is an arbitrary function of a scalar field on the brane. We present 
four-dimensional equations on a DGP brane with a scalar field non-minimally coupled 
to the induced Ricci curvature, embedded in a five-dimensional Minkowski bulk. This 
is an extension to a braneworld context of scalar-tensor (Brans-Dicke) gravity. We show 
that our model allows for an embedding of the standard Friedmann cosmology in the 
sense that the cosmological evolution of the background metric on the brane can entirely 
be described by the standard Friedmann equation plus energy non-conservation on the 
brane. In this framework we explore the relation between our formalism and the so-called 
dark radiation formalism. We study cosmological implications of both minimal and non- 
minimal extension of our model. In minimal case, motivated by modified theories of 
gravity, the potential describing minimally coupled scalar field is taken to be that with 
R^ 1 term added to the usual Einstein-Hilbert action[ll]. In non-minimal case however, 
we concentrate mainly on the potential of the type V((j>) = \(f) n . We study the weak field 
limit of our model and show that the mass density of ordinary matter on the brane should 
be modified by the addition of the effective mass density attributed to the non-minimally 
coupled scalar field on the brane. Considering the case of FRW brane, we obtain the 
evolution of the metric and scalar field by solving the field equations in the limit of small 
curvature. Our solutions for minimal case predict a power-law acceleration on the brane 
supporting observed late-time acceleration. For non-minimal case (by adapting a simple 
ansatz) we show that by a suitable choice of non-minimal coupling and scalar field poten- 
tial one can achieve accelerated expansion in some special cases. We study quintessence 
model with non-minimally coupled scalar field on the brane and discuss some observa- 
tional constraints imposed on the value of the non-minimal coupling using supernova data. 
We use a prime for differentiation with respect to fifth coordinate except for two special 
cases: a' = and V = An overdot denotes differentiation with respect to the 
comoving time, t. 
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2 Induced Gravity with Non- Minimally Coupled Brane- 
Scalar Field 



The action of the DGP scenario in the presence of a non-minimally coupled scalar field 
on the brane can be written as follows 



-/ 



y=0 
(1) 

where we have included a general non- minimal coupling a(<f>) in the brane part of the 
action( for an interesting discussion on the possible schemes to incorporate NMC in the 
formulation of scalar-tensor gravity see [22,35]). y is coordinate of the fifth dimension 
and we assume brane is located at y — 0. qab is five dimensional bulk metric with Ricci 
scalar 1Z, while is induced metric on the brane with induced Ricci scalar R. qab and 
q^y are related via q^ u = 5 tJ A 5 v L ^ Qab- K is trace of the mean extrinsic curvature of the 
brane defined as 



K 



- lim ( 

2 e^0\ 



K, 



+ 



K, 



(2) 



and corresponding term in the action is York- Gibbons- Hawking term[32] (see also [20]). 
The ordinary matter part of the action is shown by Lagrangian C m = £ m (g^ , ip) where 
if) is matter field and corresponding energy-momentum tensor is 



T = 



"2t : h q^i/i-m- 



The pure scalar field Lagrangian, C 4 
energy-momentum tensor 



-\q^^u<P - V((j>), yields the following 



1 



= v M 0v^0 - -<v(V0) 2 - q lu y{4>) 



(4) 



The Bulk-brane Einstein's equations calculated from action (1) are given by 

m\ (u A B ~ ^9abK) + 



ml8 A ^ B v 



a(0) ( - ) - V M V„a(0) + g^n {4) a(0) 5(y) = WV(3/), (5) 



where is 4-dimensional(brane) d'Alembertian and Y M „ = T^ + r^ . This relation can 
be rewritten as follows 



m\ ( Kab 



g AB n) + m\a{<i>)S A > t 8B V [r^ - \q^R) %) = <5a (6) 
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where T^ v is total energy-momentum on the brane defined as follows 



From (6) we find 



and 



= m\V M V ' v a(4>) - m^g^D (4) a(0) + T^, 



Gab = TIab — -9abT^ — 



Gfj,u — I R 



L fJ,U 



T 

1 \w 



mla(<f>) 



(7) 

(8) 
(9) 



for bulk and brane respectively. The corresponding junction conditions relating the ex- 
trinsic curvature to the energy-momentum tensor of the brane, have the following form 



lim 





y=+e 1 








v=-e ml . 



m\a{(j>) 



y=0 



ml 



Ruu ~ 7:QnuQ a>3 Ra/3 

D 



(10) 



We set qab = Vab + ^ab where Hab are small perturbations to investigate the 
weak field limit of the scenario. Within the Gaussian normal coordinate, we impose the 
following harmonic gauge on the longitudinal coordinates [9] 



ft h a + ft h — -ft (h a + h ) 



(11) 



This will led us to a decoupled equation for the gravitational potential of a static mass 
distribution. The transverse equations in the gauge (11) give h yjX = 0, h yy = h a a . The 
remaining equations take the following form 



m 3 4 (d a d a + d 2 y ) V + m\a{<f>) (d^h^ - d^h a a ) 5{y) = -25{y) 



(12) 

We suppose that non-minimally coupled scalar field has an effective mass M^. The grav- 
itational potential of mass densities p$(r) = M$6(r) and p${r) = M ( / ) 5(r) on the brane 
satisfies the following equation 



m, 



(d a d a + &l) U(f,y) + mla(cf))8(y)d a d a U(r,y) = -{M* + M^)5(r)5(y). (13) 



This equation shows that in the presence of non-minimally coupled scalar field, the mass 
in standard DGP framework should be modified by the addition of the mass of the non- 
minimally coupled scalar field. Using the following Fourier ansatz 



U(f,y) = -^-^ J d 3 p J dp y U(p,p y )exp [i{p-f+p y y) 



(14) 
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in equation (13), we find 



m\{f + p 2 y )U(p,p y ) + H^-f J d P ' y U(p,p' y ) = - 2 -{M,, + M+). (15) 



This integral equation has the following solution 

4 + M 



U(j>,Py) = —. 



3 (P 2 + P 2 y) i^ m l + m 3«(0) \P\) 

The resulting potential on the brane is 



(16) 



2 2 
cos(£ a r) cos(£ Q r) Si(^r) + - sin(£ a r) Ci(£ a r) 

71 71 



(17) 



where £ Q 



2m| 



2 and the sine and cosine integrals are defined as follows 

r°° cos a; 



Si(x) = / 



21 , smw 
aa; . 

UJ 



Ci(x) — — I duo 

Jx to 



Now there is a modified transition scale 



^1«(0) 

2m? 



a(0)£ 



DGP 



(18) 



between four and five- dimensional behavior of the gravitational potential in this scenario: 



r « £ a : f/(f) = 



and 



r > 4 : 



6/T77l3Q;(0)r 



C/(r) 



2 \ r / t r ^ 



7lH a 



§7x 2 m\r 2 



(19) 



(20) 



where 7 = 0.577 is Euler's constant. Therefore, existence of non-minimally coupled 
scalar field on the brane has the following consequences: the mass density of ordinary 
matter on the brane should be modified by the addition of the mass density attributed 
to the scalar field on the brane and the DGP transition scale between four and five 
dimensional behavior of gravitational potential now is explicitly dependent on the strength 
of non-minimal coupling. If a(<f)) varies slightly from point to point on the brane, it can 
be interpreted as a spacetime dependent Newton's constant. The dynamics that control 
this variation are determined by the following equation 



, dV ml /da((f)) 



)r = 0. 



(21) 
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Note that in the real world we don't want a(<f>) to vary too much since it will have ob- 
servable effects in classic experimental tests of general relativity and also in cosmological 
tests such as primordial nucleosynthesis [33]. This can be ensured either by choosing a 
large mass for scalar field 4> or choosing a{4>) so that large changes in give rise to rel- 
atively small changes in Newton's constant. So, when a(4>) varies in DGP brane from 
point to point, the crossover scale will change and is no longer a constant. This feature 
would change previous picture of crossover scale in DGP scenario and my change some 
arguments on phenomenology of this scenario [34] . Note that we can define a modified 
four-dimensional Planck mass as = \J a{(p)m^ and therefore the effect of non-minimal 
coupling can be attributed to the modification of four dimensional Planck mass. This is 
equivalent to modification of four dimensional Newton's constant [22]. If we use the re- 
duced Planck mass for 777,3, gravitational potential for small r limit will differ from the 
ordinary four- dimensional potential by a factor In fact the coupling of the masses 

on the brane to the induced Ricci tensor on the brane is modified by this factor. Existence 
of this extra factor is in agreement with the tensorial structure of the graviton propagator 
due to additional helicity state of the five- dimensional graviton [9,34]. Figure 1 shows 
the shape of the DGP potential for r <ti £ a (equation (19)) with some arbitrary values 
of the NMC. As this figure shows, for large negative values of NMC, it is possible to 
have repulsive gravitational potential. For positive values of NMC the DGP potential for 
r <C i a is always attractive. Figure 2 shows the DGP potential for r ^> £ a (equation (20)). 
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Figure 1: DGP potential for r <^ £ a and with some arbitrary values of NMC (equation (19)). 
Ignoring the tensor structure, the gravitational potential of a source of mass m is given 
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Figure 2: DGP potential for r » l a and with some arbitrary values of NMC (equation (20)). 

by Ugrav ~ - G branem f or r <^ £ a an d U grav ~ -QmmSl f or r > £ a respectively. That is, 
the potential exhibits four- dimensional behavior at short distances and five-dimensional 
behavior (i.e., as if the brane were not there at all) at large distances. In the absence of 
NMC, for the crossover scale to be large, we need a substantial mismatch between four- 
dimensional Planck scale (corresponding to the usual Newton's constant, G\, rane = G) and 
the fundamental, or bulk, Planck scale M 4 [34]. However, with NMC a((f>), evidently £ a 
can be large with large values of a((f>). Therefore existence of NMC changes the mismatch 
condition between four- dimensional Planck scale and the fundamental Planck scale. In 
minimal case, the fundamental Planck scale M 4 has to be quite small in order for the 
energy of gravity fluctuations to be substantially smaller in the bulk versus on the brane, 
the energy of the latter being controlled by M 3 = M P . This situation can be mediate due 
to the presence of NMC. Note that when M 4 is small, the corresponding Newton's constant 
in the bulk, Gb u ik, is large. So, for a given source mass m, gravity is much stronger in the 
bulk. Earlier work using Supernova data implies that the best-fit for density parameter of 
ordinary matter is given by tt° M = 0-l8t^l [41,42]. We may introduce a new effective dark 
energy component, Qe a , where Qe a = to resort the identity: 1 = flM + ^£ a - Therefore 
we find £ a = (l. 2 l+o^g) if J" 1 . Assuming a flat universe, a more recent supernova data [43] 
suggests a best-fit of Vl M = 0.21 corresponding to a best-fit crossover scale £ a = 1.26/fJ" 1 . 
Since 4 = or a(0) = H^, for M 3 ~ 10 18 GeV and M 4 ~ 10 3 GeV, we find 

a(4>) ~ 2.52 x IQ-^Hq 1 . 
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3 Non- Minimal DGP Cosmology 



It has been shown that DGP model can account for the standard Friedmann cosmology 
at any distance scale on the brane [20]. In our proposed framework, we start with the 
following line element to derive cosmological implications of our model, 

ds 2 = q^dafdx" + b 2 (y, t)dy 2 = -n 2 (y, t)dt 2 + a 2 (y, t)-f ij dx i dx j + b 2 (y, t)dy 2 . (22) 



In this relation 7^ is a maximally symmetric 3-dimensional metric defined as 



7jj ^ij ^~ 



X i X j 



kr 2 



(23) 



where k = —1,0, 1 parameterizes the spatial curvature and r 2 = XiX 1 . We assume that 
scalar field (f> depends only on the proper cosmic time of the brane. Choosing gauge 
b 2 (y,t) = 1 in Gaussian normal coordinates, the field equations in the bulk are given by 
(8) with the following Einstein's tensor components 



_ 2 , a 2 a' 2 a" k \ 

Lt — ^ n [ o 9 o" 1 9 J ) 



Gij = 'JijO^ 



/2 

/a 



a 2 
n 2 a 2 



n A ' \ n 



a" n'a' 



a na n 
a ' na n 2 a n 3 a 2n 



G 



Oy 



n' a a' 



n a a 



G 



yy 



a 2 k n'a' na d 
n 2 a 2 a 2 na n 3 a n 2 a 



3(V 



The field equations on the brane are given by the following equations 

^ fc\ 2 



Gffi = 3n 2 



na 



a 



+ 



n 3 a n 2 a 



Gij — 7ij 
and scalar field evolution equation 

, „d n 

(/>■ 



n 2 a 2 a 2 J m§a(0) 



k 



00 • 



+ n 2 



n 2 a 2 ~^ a 2 



,dV m 



T 



a n, 

where Ricci scalar on the brane is given by 



- '-fn'a'Riq] = 0, 



a 2 n 2 



6^ + 6H 2 -6^ 
a K a J a 



.an 

, n 



(24) 

(25) 
(26) 
(27) 

(28) 
(29) 

(30) 
(31) 
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The other important equation is the continuity equation on the brane. Suppose that 
ordinary matter on the brane has an ideal fluid form, T^ u — (p + p)u^u v + pq^. Since 
K tt = nn' and K rr = —ad ', equation (10) gives the following matching conditions 



m 



m| 



a{4>) 



n 2 a 



+ 



k 



lim n [^n]^:(t) = -|(2n) 



2 

mi 



j/=0 
2 



(p + p<^)« 

3m? 



(32) 



y=0 



/ a a na k \ 



+ 



j/=0 



3m? 



2(p + ^) + 3(p + ^) 



(33) 



j,=0 



where energy density and pressure of non-minimally coupled scalar field are given as 
follows 

6a'H(j> 



\^ + n 2 \/(0) 



H 

n z 



2a" 



(34) 
(35) 



J y=o 



and if = - is Hubble parameter. Note that part of the effect of non-minimal coupling of 
the field is hidden in the definition of the effective energy density and pressure which 
both include non- minimal terms. Now using (26), since in the bulk G 00 = 0, we find 



lim 



~n'~ 


y=+e 


a'- 


y=+e 


. n _ 


y=-e 


. a _ 


y=-e 



(36) 



using relations (32) and (33) we find the following relation for conservation of energy on 
the brane 

p + P4, + 3#(p + P4> + p + p^) = 6a>(# 2 + — ) . (37) 

a 

Thus the non-minimal coupling of the scalar field to the Ricci curvature on the brane 
through a(4>) leads to the non-conservation of the effective energy density. 

To obtain the cosmological dynamics, we set n(0,t) = I. With this gauge condition 
we recover usual time on the brane via transformation t = /*n(0, r))dr). In this situation, 
our basic dynamical variable is only a(y, t) since n(y, t) now is given by 

a(y,t) 



n(y,t) 



a(0,t)' 



(3? 



where H = ^jjj^y is Hubble parameter on the brane. Now we can write the basic set of 
cosmological equations for a FRW brane in the presence of a non-minimally coupled scalar 
field. The first of these equations is given by matching condition 



(*) = -§ 
m% 



a(4>) 



n 2 a 



k 
a 



y=0 



(P + P4>) a 
3mi 



(39) 



y=0 
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Insertion of ^ = 4- into equations (24) and (27) yields 

^-a'a 3 Goo = ^r{%a 2 ~ a'V + ka 2 
6n A ay \n z 

and 



• 3 ft 
— (10, (jyy 



d fa 2 



2 II 2 , 7 2 

a — a a + /ea 



<9t \n 2 

These two equations imply that the Binetruy et al [14] integrals 



0. 



x + 



and 



T~ = 



j/>0 



j/<0 



(40) 



(41) 



are constant and if a' is continuous on the brane then X + = X . These equations along 
with scalar field equation 



, a n\ ■ 9 dV 9 da , 

+ 3 U + n 2 — - n 2 —R[q] = 0, 

an/ ' 



and 



n(y,t) 



d(f> d(j> 
a(0,t)' 



(42) 



(43) 



constitute the basic dynamical equations of our model. In the absence of transverse 
momentum, T Qy = 0, one can show that X + = X~. In fact I ± can be considered as 
initial conditions and these quantities reflect the symmetry across the brane. In the case 
of X + 7^ X~ there can not be any symmetry across the brane. So we first consider the 
case X + = X~ in which follows. Our cosmological equations on the brane now take the 
following forms(note that n(0,t) = 1) 

d 2 (0,t)+k (p + P<p ) 



a 2 (0,t) 3m|a(0) : 



+ 3 
X = 



a(0,t) i dV(<f>) 



a(0, t) d(f) 
d 2 (0,t)-a' 2 (y,t) + k]a 2 (y,t) 
a(y,t) 



n(y,t) 



a(0,t)' 



Using equation (46), the scale factor is calculated as follows 



a 2 (y, t) = a 2 (0, f) + d 2 (0, t) + k y 2 + 2 



(d 2 (0,t) + k)a 2 (0,t) -X 



(45) 
(46) 
(47) 

(48) 
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and therefore n(y,t) is given by equation (47); 



n(y, t) = a(0, t) + a(0, t)y 2 + a(0, t) 



a(0,f)a(0,t) +d 2 (0,t) + k 
(a 2 (0,t) + k)a 2 (0,t) - I 



x 



a 2 (0,t) + [a 2 (0,t) + k\y 2 + 2[(a 2 (0,t) + k)a 2 (0,t) - I 



(49) 



So, the components of 5-dimensional metric (22) are determined. If we set initial condi- 
tions in such a way that X = 0, we find the following simple equations for cosmological 
dynamics 



a(y,t) = a(0,t) + \a 2 (0,t) + k 
d(0,t) 



n(y,t) = 1 + 



--y- 



(50) 
(51) 



^Ja 2 (0, t) + k 

Therefore, our model allows for an embedding of the standard Friedmann cosmology in 
the sense that the cosmological evolution of the background metric on the brane can be 
described by the standard Friedmann equation plus energy non-conservation on the brane. 

So far we have discussed the case X + = X~ with a continuous warp factor across the 
brane. In the case of X + ^ X~, there cannot be any symmetry across the brane. In 
this case the basic set of dynamical equations is provided by equations (39), (40) and 
(41) plus the non-conservation of the effective energy density given by (37). In this case, 
evolution of the scale factor on the brane is given by elimination of a'(y — > ±0,t) from 
the following generalized Friedmann equation 



a 2 (0,t) + k-a- 2 (0,t)l + 



/ ^m 2 f a 2 {0,t)^ 



a 2 {0,t) + k-a- 2 (0,t)l~ 
k\ (p + P4,)a(0,t) 



3m? 



(52) 



This is the most general form of the modified Friedmann equation for our non-minimal 
framework. After determination of a(0, t), since X ± are constants, a(y,t) can be calcu- 
lated from (46). This is the full dynamics of the system. Note that in the case where 
the right hand side of equation (52) is negative, at least one sign in the left hand side 
should be negative depending on initial conditions. However, the dynamics of the problem 
does not require symmetry across the brane. Therefore, we have shown the possibility of 
embedding of FRW cosmology in DGP scenario with a 4D non- minimally coupled scalar 
field on the brane and equation (52) is the most general form of FRW equation in this 
embedding. 
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4 Dark Radiation Formalism with Non-Minimal Cou- 
pling 

In order to show that our model is consistent with dark energy formulation of Friedmann 
equation, we first obtain non-minimal extension of this formalism and then we show 
that this equation can be obtained in our framework more easily. The effective Einstein 
equation on the brane is given by [24] 

G^ u = g &yuvi (53) 



where 



n„ = --J^Tf + iTT^ + \g» v (T pa To° - It 2 ), (54) 



and T^ v is given by equation (7). Also we have 

S,u = C MRNS n M n N g R ^ g s u (55) 

where Cmrns is fi ye dimensional Weyl tensor and ua is the spacelike unit vector normal 
to the brane. Now, using equation (53) we find 

G% = ^ - £°o (56) 



where for FRW universe we have 



G% = -3(H 2 + ^). (57) 
Similarly, for space components we have 

IT 

where 

G l , = -(2ij + 3# 2 + ^)<5Y (59) 

Now using equation (54) we find n° = -^( T °o) 2 and = -^°o(^°o - 2T 1 1 )S i j . 
Also, equation (7) gives 

T°o = -(p + p^)-m 2 a(0)G' O o (60) 

and 

T'^-ip + p^i-mla^G'j, (61) 
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where p<j, and are given by (34) and (35) with n(0,t) = 1. These equations lead us to 
the following generalized Friedmann equation 



3(h 2 + ^)=£\ 



+ 



p + Pt-3mla(<j)){H 2 + —) 

a 



(62) 



Using Codazzi equation we have V"^ = (see for example [24,31]). Therefore we find 
S°q + 4:HS°o = which integration gives S° = ^ with £ as an integration constant. 
Therefore, equation (62) can be re-written as follows 



H 2 + — 



a 2 3m 2 a(<f)) 









(p+p 


<p + Po 


l + eJl + — 


L V Po l 



p + p<p -m§a(0)-j 



(63) 



where po 



m|a(</>) 



and e = ±1 shows the possibility of existence of two different 



branches of FRW equation. In the high energy regime where ^> 1, we find 



H 2 



(p + P<t> + e \j 2 (p + P4>)Po 



3m 2 a(4>) 



(64) 



which describes a four dimensional gravity with a small correction. In the low energy 



regime where < 1, we find 



Po 



H 2 « 



(1 + e)(p + P(j) ) + (l + e) po 



g (P + P<^ 
4 Po 



(65) 



Now we show that this result can be obtained in our formalism more easily. For this 
purpose we show that relation (52) for the case with X + = Z~ = X and a discontinuous 
warp factor across the Z 2 symmetric brane leads to this result with some simple algebra. 
For simplicity, we define x = H 2 + b = p + p^, y = a(0)m|, and z = m^. With 
these definitions, equation (52) (with upper sign for instance), transforms to the following 
form 



J 4 



x — 



+ U - 



-X- — 

z 3z 



(66) 



Solving this equation for x (with I + = 1 = J) gives the following result 



x = 



(67) 



A little algebraic manipulation gives 

1 



x = 



3y 



6z 2 6z 2 
6H ± 

y y 



'i + 



by ly 2 



3z 2 a*z 2 



(68) 
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Considering both plus and minus signs in equation (52) and using original quantities we 
obtain 



P + P<t> + Po 



2 S 
l + eJl + — p + p^-m§a(0)-J 
V Po 1 a 



(69) 



where p = = m^ a \^) an d &o = 3Z is a constant. This analysis shows the consistency 
of our formalism with dark-radiation formalism presented above. 



5 Cosmological Considerations 

Now to discover cosmological implications of non-minimally coupled scalar field on DGP 
braneworld, we proceed as follows: the evolution of scalar field on the brane for spatially 
flat FRW metric (k = 0) is given by the following equations 

*° = S(kl (p + ^ (70) 

and 

4 + 3H i+?^ = a'R[q], (71) 
d(p 

where H = and for n(0,t) = 1 (on the brane) we have 

P4> = \tf + V{4>) - Qa'H<j). (72) 
Ricci scalar on the brane is given by the following relation 

R = 6- + 6[-) . 73 

n . v nJ 



In the absence of ordinary matter on the brane we set p = 0. To find cosmological 
implications of our model we should solve equations (70) and (71) using (72) and (73). 
For comparison purpose, we first give an overview of the minimal case, that is, the case 
with a = constant which has been discussed with details in reference [30]. First we 
should specify the form of the scalar field potential V(4>). We use the following potential 
which has been motivated from theories of modified gravity with Lagrangian of the type 
C(R) = R- £ with R" 1 term[ll,30] 

V^^^mle^l-^). (74) 
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In this case we find a(0, t) oc t 4 / 3 and oc — | \nt for the evolution of scale factor and 
scalar field on the brane respectively. Obviously, this predicts a power-law acceleration on 
the brane. This result is consistent with the observational results similar to quintessence 
with the equation of state parameter —1 < wde < — §, [6,20]. In this case, using (48) 
and (49), the evolution of a(y,t) and n(y,t) is given by the following equations [30] 

16 2 \ /16 



a\y,t) = C 2 (tl + ftly 2 ) + 2{^CHf - x) % (75) 



and 



n(y,t) = C 
where C is a constant and X 



-i 



1 4._2 2 20 2 8 /16 4 in _\ 7 



' (76) 



a(y,0' 

2 . If we set the initial conditions 



d 2 (0,t) - a ,2 (y,t) + k 
in such a way that X = 0, these results become very simple 

a(y, t) = C(tl + U>y) , n(y, t) = C(l + |) . (77) 

In the case where a(4>) is not just a constant, the situation becomes more complicated 
since now the kind of the cosmological solutions depend explicitly on the value of non- 
minimal coupling. We first try to obtain a necessary condition for the acceleration of the 
universe in an specific model. Suppose that p = 0. By definition, the required condition 
for acceleration of the universe is p^ + 3^0 < 0. Using equations defining p^ and p^, we 
find 

(1 + 3a")0 2 - V(<f>) + 3a'(H<p + 0) < 0. (78) 
Using Klein-Gordon equation (71), this relation can be rewritten as follows 

(1 + 3a")0 2 - VU) + 3a' 2 R - 6a' H <j) - 3a'^- < 0. (79) 

Finally, using (72), this can be written as 

P4> - 2V(<t>) + (- + Sa")<p 2 + 3a' 2 R - 3a'^- < 0. (80) 
2 deb 

This is a general condition to have an accelerating universe with non-minimally coupled 
scalar field on the brane [22,35]. To proceed further, we assume weak energy condition 
p^ > 0. Motivated from several theoretical evidences( for example: conformal coupling 
in general relativity and other metric theories [35], renormalization group study of A0 4 
theory[36] and large N limit of the Nambu-Jona-Lasinio model[37], (see also [22] and 
references therein), in which follows we set cx(<f>) = |(1 — £0 2 ) with £ < |. Therefore we 
find 

3£ dV 

V - > 0. (81) 
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As an example, suppose that V(<f>) = \<p n . In this case with £ < | and using (81), we 
find 

(82) 



A(l-^]>0. 



If we assume a positive scalar field potential with A > 0, the condition for accelerating 
expansion restricts £ to the values which £ < So in the presence of non- minimally 
coupled scalar field, it is harder to achieve accelerating universe with usual potentials. 

Now with above definitions of non-minimal coupling and scalar field potential, equa- 
tions (70) and (71) take the following forms 

a 2 2 



■(l-£0 2 ) 



2\-l 



U 2 + xr + 6^-0 

2 a . 



and 



+ 3-0 + n\(t) n - v + 6£0 
a 



a a 

- + - 

a o, 



0. 



(83) 
(84) 



In order to study late time behavior of these equations, we try the following ansatz, 

a(t)«Af, Br" (85) 

where we have assumed a decreasing power law ansatz for scalar field. With these choices 
and setting n — 4, equation (83) gives 



3z/_ 2 

2 1 2 -1 



- -^Vr 2 "" 2 = (-B 2 fi 2 - Q^B 2 fiu)r 2 ^ 2 + A£ 4 r 4 " 
On the other hand, equation (84) gives 

fi(fi + 1) - 3fiu + 6£(2z/ 2 + v)]t-»- 2 + 4A5 3 r 3 ^ = 0. 



(86) 



(87) 



Considering terms of order (9(t^~ 2 ), equation (86) gives £ > ^ if we require a positive 
and real v. So for this special ansatz £ is restricted to the condition ^ < £ < |. With 
the same procedure, equation (87) gives 

/i 2 + (1 - 3z/)/i + 12£z/ 2 + 6£z/ = 0. (88) 

Again, positivity and reality of solutions for \x lead to the following constraint 

(9-48£)^ 2 - (6 + 24£)^+ 1 > (89) 

where for £ = ^ and taking equality, we find v = 4± ~f^ which plus sign obviously leads 
to an power-law accelerated expansion. So, we conclude that although in the presence of 
non-minimally coupled scalar field, accelerated expansion of universe is harder to achieve 
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-02 -0.1 0.1 0.2 

e=7 

£=o 

£=-.7 

Figure 3: An effective potential V"(0) in the Einstein frame in the case of £ = 0.7, 0, —0.7 with 
p = 1/2 (top, middle, bottom). When £ is positive, since the potential V(4>) is flat in the region 
of 4> > compared with the £ = case, we can expect assisted inflation to occur in this region. 
When £ is negative, assisted inflation can be realized in the region of <j) < 0. 



relative to minimally coupled scalar field case, but with a suitable choice of nonminimal 
coupling it is possible to explain this accelerated expansion. 

On the other hand, with a nonminimal coupling of the kind — ^R(p 2 , we can define 

Geff = (90) 

2 

where 2 = In order to connect to our present universe, G e ff needs to be positive 
for the case of the positive £, which yields \4>\ < 4> c = j^= . When £ is negative, such a 

constraint is absent. By performing a conformal transformation [40] to transform to the 
Einstein frame, we define 

v s m ox) 
- (i-^0 2 ) 2 

Figure 3 shows the behavior of this effective potential with different values of nonminimal 
coupling and V(<f)) defined in (74). As we see, when £ is positive, since the potential 
V(4>) is flat in the region of (ft > compared with the £ = case, we can expect assisted 
inflation (due to NMC) to occur in this region. When £ is negative, assisted inflation can 
be realized in the region of < 0. 

Based on some holographic dark energy model, the value of NMC, £ to have an ac- 
celerated universe is restricted to the interval 0.146 ^ £ ^ 0.167 [44]. On the other hand, 
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current experimental limits on the time variation of G constraint the nonminimal cou- 
pling as — 1CT 2 ^ £ r> 10~ 2 [45]. Solar system experiments such as Shapiro time delay 
and deflection of light [33] have constraint Brans-Dike parameter to be ujbd > 500. This 
constraint leads to the result of |£| ^ 2.2 x 10™ 2 for non-minimal coupling [45]. 



6 Quintessence Model 

To study quintessence model within our setting, we write the equation of state for scalar 
field which takes the following form 

V(j> 2 -2V(4>) +4(a'4) + 2Ha'<j) + a"<j) 2 ) 
W ~ P4> ~ 2 + 2V(<j>) - \2a'H4> ' ^ 

Causality implies that \w\ < 1. When = 0, we obtain p^ = —p$. In this case p$ is 
independent of a and V(4>) plays the role of a cosmological constant. If V{4>) = 0, we find 

2 + 4(a'0 + 2Ha'6 + a"dA 

w = ^ ■. '-. (93) 

2 - 12a'H(f) 

In the case of minimal coupling (a = const.) this is corresponding to a massless scalar 
field which plays the role of a stiff matter since p^ ~ ^. However, in our case with non- 
minimally coupled scalar field, the situation is very different since now a plays a crucial 
role and depending on the form of a we may obtain some traces of stiff matter or fail to 
have such extreme case. So, in principle our model provide a mechanism to avoid stiff 
matter. In the minimal case when 2 < V((f>), we obtain < — ^ which shows a late- 
time accelerating universe. In this case, since V(<f>) = \{l — w)p ( j > and = + w) In a, 
the following potential which is a Liouville-type potential, decreases when scalar field <f> 
increases and therefore gives the required quintessence 

V(<f>) = V exp ( - ^3(1 + ^)0) , (94) 

where 0<3(l + w;)<2. On the other hand, for non- minimal case we find 

V(4>) = -a(l - w)p + (a'0 + 5Ha'<p + a"0 2 ) . (95) 

In this case the quintessence potential has an explicit dependence on the non-minimal 
coupling and its derivatives with respect to 0. To have a quintessence model, we should 
impose some limits on the shape of this non-minimal coupling. Since based on quintessence 
proposal the dark energy of the universe is dominated by the potential of scalar field, 
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we should impose suitable constraints on non-minimal coupling such that scalar field 
potential decrease when scalar field increases. We assume ordinary matter on the brane 
has energy density p oc a~ 3 and vanishing pressure p — 0. Now total energy density 
becomes pr = p<j> + p- The necessary and sufficient condition for acceleration of the 
universe, p T + 3P T < 0, leads to the following relation 

(1 + 3a")0 2 - V{4>) + 3a' 2 R - 6a'# o - 3a'^f- + P - < 0, (96) 

dcp 2 

which can be written as 



p - 2V{4>) + Q + 3a")0 2 + 3a' 2 R - 3a'^- + P - < 0. (97) 

Z (X(p z 

This is a general constraint to have a quintessence scenario in the presence of non-minimal 
coupling. If we assume that p and p^ are non- negative, with a(<f>) = \{l — £0 2 ) and £ < |, 
we find 

V-^->0. (98) 
To have quintessential expansion, this constraint with < £ < \, restrict the form of 



scalar field potential to potentials V(<j>) > where ^ 

oo = ^(l-J^\t=tn), = p/Pr, fiwi = PrA/Pr. and £ 

14, 



In 



ft(f)"-p(-f 



< where 



3?v 2n i*=*oy' 



VI = p/ p c , £l<t> = P<p/Pc and has been approximated by its present 
value [22]. So, to have a quintessential expansion with non-minimal coupling with < £ < 
1/6 we need a potential that does not grow faster than f((f>) = Vo(^) exp (^-) with 
variation of 0. An inverse-power-law potential such as V(<f>) = p 4+6 (fi~ s , an exponential 
potential such as V(4>) = p 4 exp(— A^-) and several other possibilities[38,39] provide 
suitable potentials for quintessence. But obviously quintessence model in the presence of 
non-minimal coupling needs more artificial arguments than minimal case. 



7 Summary and Conclusions 

In this paper we have considered the DGP model with a non-minimally coupled scalar 
field on the brane. As we have explained, the introduction of non-minimal coupling is not 
just a matter of taste; it is forced upon us in many situations of physical and cosmological 
interests such as quantum corrections to the scalar field theory and its renormalizability 
in curved spacetime. In the spirit of DGP inspired gravity, we have studied the effect of an 
induced gravity term which is an arbitrary function of a scalar field on the brane. We have 
presented four-dimensional equations on a DGP brane with a scalar field non-minimally 
coupled to the induced Ricci curvature, embedded in a five- dimensional Minkowski bulk. 
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This is an extension to a braneworld context of scalar-tensor (Brans-Dicke) gravity. Cos- 
mological implications of both minimal and non-minimal extension of our model are stud- 
ied. In minimal case, we have considered an exponentially decreasing potential which has 
been motivated by modified theory of gravity with R' 1 modification. In the non-minimal 
case however, we have considered potential of the type V{4>) = \<p n . We have stud- 
ied the weak field limit of our model and it has been shown that the mass density of 
ordinary matter on the brane should be modified by the addition of the effective mass 
density attributed to the non-minimally coupled scalar field on the brane. Also in this 
case crossover scale of DGP scenario is modified by the presence of non-minimal coupling. 
We have discussed the role of nonminimal coupling in crossover distance. Considering 
the case of FRW brane, we have obtained the evolution of the metric and scalar field by 
solving the field equations in the limit of small curvature. Our solutions for minimal case 
predict a power-law acceleration on the brane supporting observed late-time acceleration. 
For non-minimal case we have shown that by a suitable choice of non-minimal coupling 
and scalar field potential one can achieve accelerated expansion in some special cases. 
However, As Faraoni has shown, in the presence of non-minimally coupled scalar field 
accelerated expansion of universe is harder to achieve relative to minimally coupled scalar 
field case. We have studied quintessence model in our framework and it has been shown 
that for a restrict class of non-minimal coupling one can achieve quintessence potential. 
As an important achievement our analysis shows that DGP model allows for an embed- 
ding of the standard Friedmann cosmology in the sense that the cosmological evolution 
of the background metric on the brane can entirely be described by the standard Fried- 
mann equation plus energy non-conservation on the brane. As we have shown, our model 
gives dark energy extension of Friedmann equation more easily than standard framework. 
Some observational and experimental constraint on nonminimal coupling are discussed. 
The issue of Non-minimal inflation on warped DGP braneworld and confrontation with 
WMAP3 data will be reported in a separated paper [46]. 

Acknowledgement 

It is a pleasure to appreciate members of the Centre for Particle Theory at Durham 
University, specially Professor Ruth Gregory for hospitality during my sabbatical leave. 
Also I would like to appreciate referee for his/her important contribution in this work. 



22 



References 



[1] A. G. Riess et. al. (Supernova Search Team Collaboration), Astron. J. 116 (1998) 
1006, [astro-ph/9805201]. 

[2] S. Perlmutter et. al, Astrophys. J. 517 (1999) 565, [astro-ph/9812133], 

D. N. Spergel et. al, Astrophys. J. Suppl. 148 (2003) 175, [astro-ph/0302209]. 

[3] C. L. Bennett et. al, Astrophys. J. Suppl. 148 (2003) 1, [astro-ph/0302207]. 

[4] C. B. Netterfield et. al, Astrophys. J. 571 (2002) 604, [astro-ph/0 104460], 
N. W. Halverson et. al, Astrophys. J. 568 (2002) 38, [astro-ph/0 104489]. 

[5] D. N. Spergel et al, WMAP Three Year Results: Implication for Cosmology, 
arXiv:astro-ph/0603499 

[6] T. Padmanabhan, Phys. Rept. 380 (2003) 253 

V. Shani and A. A. Astarobinsky, Int. J. Mod. Phys. D 9 (2000) 373 
S. M. Carroll, Living Rev. Rel. 4 (2001) 1 

E. J. Copeland, M. Sami and S. Tsujikawa, arXiv:hep-th/0603057 
S. Nesseris and L. Perivolaropoulos, arXiv:astro-ph/0610092 

[7] R. R. Caldwell, Phys. Lett. B 545 (2002) 23 

R. R. Caldwell, M. Kamionkowski and N. N. Weinberg, Phys Rev. Lett. 91 (2003) 
071301 

J. M. Cline, S. Y. Jeon and G. D. Moore, Phys. Rev. D 70 (2004) 043543 

[8] R. R. Caldwell, R. Dave and P. J. Steinhardt, Phys. Rev. Lett. 80 (1998) 1582 
P. J. E. Peebles and A. Vilenkin, Phys. Rev. D 59 (1999) 063505 
P. J. Steinhardt, L. M. Wang and I Zlatev, Phys. Rev. D 59 (1999)123504 
M. Doran and J. Jaechel, Phys. Rev. D 66 (2002) 043519 
A. R. Liddle, P. Parson and J. D. Barrow, Phys. Rev. D 50 (1994) 7222 

[9] G. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. B 485 (2000) 208, [hep- 
th/0005016] 

G. Dvali and G. Gabadadze, Phys. Rev. D 63 (2001) 065007, [hep-th/0008054], 

G. Dvali, G. Gabadadze, M. Kolanovic and F. Nitti, Phys. Rev. D 65 (2002) 024031, 

[hep-th/0106058] 



23 



[10] C. Deffayet, Phys. Lett. B 502 (2001) 199, [hep-th/0010186] 

C. Deffayet, G. R. Dvali and G. Gabadadze, Phys. Rev. D 65 (2002) 044023, [astro- 
ph/0105068] 

C. Deffayet and S. J. Landau, J. Raux, M. Zaldarriaga and P. Astier, Phys. Rev. D 
66 (2002) 024019, [astro-ph/0201164], 

J. S. Alcaniz, Phys. Rev. D 65 (2002) 123514, [astro-ph/0202492] 

D. Jain, A. Dev and J. S. Alcaniz, Phys. Rev. D 66 (2002) 083511, [astro-ph/0206224] 
A. Lue, R. Scoccimarro, G. Starrkman, Phys. Rev. D 69 (2004) 044005, [astro- 
ph/0307034]. 

[11] S. M. Carroll, V. Duvvuri, M. Trodden, M. Turner, Phys. Rev. D 70 (2004) 043528, 
[astro-ph/0306438]. 

[12] N. Arkani-Hamed, S. Dimopoulos, G. Dvali, Phys. Lett. B 429 (1998) 263, [hep- 
ph/9803315] 

N. Arkani-Hamed, S. Dimopoulos, G. Dvali, Phys. Rev. D 59 (1999) 086004, [hep- 
th/9807344] 

[13] L. Randall, R. Sundrum, Phys. Rev. Lett. 83 (1999) 3370 
L. Randall, R. Sundrum, Phys. Rev. Lett. 83 (1999) 4690 

[14] P. Binetruy, C. Deffayet, D. Langlois, Nucl. Phys. B 565 (2000) 269, [hep-th/9905012] 
P. Binetruy, C. Deffayet, U. Ellwanger and D. Langlois, Phys. Lett. B 447 (2000) 285, 
[hep-th/9905012] 

R. Maartens, Living Rev. Relativity, 7 (2004) 7, http://www.livingreviews.org/lrr- 
2004-7 

D. Wands, arXiv: gr-qc/0601078 and references therein. 

[15] I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos, G. Dvali, Phys. Lett. B 436 (1998) 
257 

D. Cremades, L. E. Ibanez and F. Marchesano, Nucl. Phys. B 643 (2002) 93, [hep- 
th/0205074] 

C. Kokorelis, Nucl. Phys. B 677 (2004) 115, [hep-th/0207234] 

[16] W. Muck, K. S. Viswanathan, I. V. Volovich, Phys. Rev. D 62 (2000) 105019, [hep- 
th/0004017]. 

[17] R. Gregory, V. A. Rubakov, S. M. Sibiryakov, Class. Quan. Grav. 17 (2000) 4437, 
[hep-th/0003109]. 



24 



[18] I. Ya. Aref 'eva, M. G. Ivanov, W. Muck, K. S. Viswanathan and I. V.Volovich, Nucl. 
Phys. B 590 (2000) 273, [hep-th/0004114]. 

[19] M. Cvetic, M. J. Duff, J.T. Liu, H. Lu, C. N. Pope and K. S. Stelle, Nucl. Phys. B 
605 (2001) 141, [peh-th/0011167]. 

[20] R. Dick, Class. Quant. Grav. 18 (2001) Rl, [hep-th/0105320] 
R. Dick, Actaphys. Polon. B 32 (2001) 3669, [hep-th/0110162]. 

[21] H. Collins and B. Holdom, Phys. Rev. D 62 (2000) 105009, [arXiv:hep-th/0003173] 
Y. V. Shtanov, arXiv:hep-th/0005193 

N. J. Kim, H. W. Lee and Y. S. myung, Phys. Lett. B 504 (2001) 323, [arXiv:hep- 
th/0101091] 

[22] V. Faraoni, Phys. Rev. D 62 (2000) 023504, [arXiv:gr-qc/0002091] 

[23] S. Mizuno, Kei-ichi Maeda and K. Yamamoto, Phys. Rev. D 67 (2003) 023516, [hep- 
th/0205292] 

[24] Kei-ichi Maeda, S. Mizuno and T. Torii, Phys. Rev. D 68 (2003) 024033, [arXiv:gr- 
qc/0303039] 

[25] S. C. Davis, JHEP 03 (2002) 058, [arXiv:hep-ph/0111351 and arXiv:hep-th/0106271] 

[26] C. Bogdanos, A. Dimitriadis and K. Tamvakis, arXiv:hep-th/0611181 

[27] A. Dimitriadis, C. Bogdanos and K. Tamvakis, Phys. Rev. D 74 (2006) 045003 

[28] K. Farakos and P. Pasipoularides, Phys. Lett. B 621 (2005) 244 

K. Farakos and P. Pasipoularides, Phys. Rev. D 73 (2006) 084012, [arXiv:hep- 
th/0610010]. 

[29] Rong-Gen Cai and H. Zhang, JCAP 0408 (2004) 017, [arXiv:hep-th/0403234] 

[30] K. Atazadeh and H. R. Sepangi, Phys. Lett. B 643 (2006) 76, [arXiv:gr-qc/0610107] 

[31] M. Bouhamdi-Lopez and D. Wands, Phys. Rev. D 71 (2005) 024010, [arXiv:hep- 
th/0408061] 

[32] J. W. York, Phys. Rev. Lett. 28 (1972) 1082 

G. W. Gibbons and S. W. Hawking, Phys. Rev. D 15 (1977) 2752 



25 



[33] S. M. Carroll, An Introduction to General Relativity: Spacetime and Geometry, Ad- 
dison Wesley, 2004 

[34] A. Lue, Physics Reports 423 (2006) 1-48, [astro-ph/0510068] 

[35] V. Faraoni, Phys. Rev. D 53 (1996) 6813 

[36] A. Bonanno, Phys. Rev. D 52 (1995) 969; A. Bonanno and D. Zappald, ibidem 55 
(1997) 6135 

[37] C. T. Hill and D. S. Salopek, Ann. Phys. (NY) 213 (1992) 21 
[38] S. Mizuno and Kei-ichi Maeda, Phys. Rev. D 64 (2001) 123521 
[39] V. Sahni, Chaos, Solitons and Fractals, 16 (2003) 527 
[40] S. Tsujikawa, Phys. Rev. D 62 (2000) 043512 

[41] C. Deffayet, G. R. Dvali, G. Gabadadze, Phys. Rev. D 65 (2002) 044023 
[42] C. Deffayet et al, Phys. Rev. D 66 (2002) 024019 

[43] [Supernova Search Team Collaboration], A.G. Riess et al, Type la Supernova Dis- 
coveries at z > 1 from the Hubble Space Telescope: evidence for past deceleration and 
constraints on dark energy evolution, Astrophys. J. 607 (2004) 665. 

[44] M. Ito, Europhys. Lett. 71 (5) (2005) 712 

[45] T. Chiba, Phys. Rev. D 60 (1999) 083508 

[46] K. Nozari and B. Fazlpour, work in progress 



26 



